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摘 要 : 准确 、 高 效 地 掌握 草原 土壤 属性 的 空间 分 
据 。 相 比 于 传统 土壤 调查 方法 ,基于 模糊 逻辑 
图 精度 ,被 广泛 应 用 于 数字 土壤 制图 
现 有 推理 模型 的 稳定 性 较 低 ,尚未 在 高 寒 草 原 
~ 地 区 约 4 km 区 域 为 研究 区 ,以 高 程 、 坡 度 、 坡 向 、 
形 因 子 为 土壤 环境 因子 ,采用 模糊 C 均值 聚 类 (下 
类 ,得 到 9 个 环境 医 
壤 一 环境 推理 方法 模拟 研究 区 表层 土壤 pH 值 空 
35 个 独立 样本 进行 精度 检验 (均匀 采样 点 16 个 ， 


草地 生态 系统 是 陆地 生态 系统 重要 的 可 再 生 自 
然 资 源 , 为 人 类 提供 气候 调节 土壤 保持 、 涵 养 水 源 
和 畜产 品 等 重要 的 服务 功能 ""。 但 目前 ,我 国 90% 
的 可 利用 草地 都 出 现 了 不 同 程 度 的 退化 ””。 已 有 
研究 表明 草地 退化 不 仅 是 地 表 植 被 的 退化 ,更 与 土 
壤 退 化 具有 密切 联系 ,植被 退化 与 草地 土壤 退化 互 
为 因果 。 因 此 , 准确 .高 效 地 掌握 草地 土壤 的 质 
量 及 空间 分 布 能 够 为 草地 资源 境 管理 提供 基础 信息 
和 参考 依据 ,对 科学 .持续 地 利用 草地 资源 具有 重要 
积极 作用 。 

传统 的 土壤 属性 空间 信息 通常 是 借助 空间 插 
值 方法 获取 ,常用 的 方法 有 普通 克 里 格 插值 法 
(Ordinary Kriging, OK ) 、 反 距离 权重 插值 法 (了 - 
verse Distance Weighted, IDW ) .普通 最 小 二 乘法 
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的 土壤 一 环境 推 到 
。 但 由 于 土壤 自身 的 空间 变异 性 及 其 与 环境 条 件 间 的 非 线 
区 进 


子 组 合 ,并 在 隶属 度 值 高 的 环境 


果 与 实测 值 基本 吻合 , 且 基 于 模糊 聚 类 和 土壤 一 环境 推理 
距离 权重 法 。 通 过 基于 模糊 逻辑 和 土壤 一 环境 推理 
证 ,结果 表明 基于 典型 点 的 采样 方案 能 够 快速 有 效 地 对 


布 能 够 为 草地 资源 境 管 理 提供 基础 信息 和 参考 依 
能够 提高 野外 采样 效率 和 预测 制 
性 ， 
高 山 草原 
数 6 个 地 
Al FR 


行 应 用 。 选 择 新 疆 巴 音 布 鲁 克 典型 3 
沿 训 面 曲 率 、 沿 等 高 线 曲 率 、 地 形 湿度 
uzzy C-means Clustering, FCM ) 方法 对 环境 
因子 组 合 中 心 共 设置 18 个 典型 点 。 运 用 土 
间 分 布 , 其 变化 范围 在 7.170 ~8.186 ZT, AH 
横 截面 采样 点 9 个 ,垂直 带 采 样 点 10 个 ) ,模拟 结 
方法 的 模拟 精度 高 于 普通 克 里 格 法 和 反 
E 的 数字 土壤 制图 方法 在 小 尺度 区 域 的 运用 验 
区 域 土壤 属性 进行 空间 模拟 ,该 方法 对 于 


已 
日 


= 类 似 小 尺度 的 研究 区 同样 有 效 。 
à 关 键 词 : 模糊 聚 类 ; 目的 性 采样 ; pH; BAAS. 


( Ordinary Least Squares , OLS) 等 ,还 有 学 者 利用 多 
元 线性 回归 分 析 集 成 所 有 成 土 因子 对 土壤 养分 进 
行 空间 分 布 预测 ,已 取得 较 多 研究 成 果 " 。 近 年 
来 , 随 着 遥感 .CS 、 人工 智能 和 模糊 推理 等 技术 的 
发 展 , 数 字 土 壤 制 图 领域 取得 了 很 大 的 进展 ”” 。 
其 中 ,基于 模糊 逻辑 的 土壤 一 环境 相似 性 推理 研 
究 日 益 受 到 关注 ,最 具 代 表 性 的 是 朱 阿 兴 等 提出 
的 SoLIM 模型 (Soil-Land Inference Model) ,其 利用 
人 工 智 能 获取 并 量化 土壤 和 环境 关系 的 知识 , 通 
过 相似 度 推理 ,使 土壤 类 型 的 推测 精度 有 较 大 提 
高 “”。 此 外 ,部 分 学 者 运用 相似 性 分 布 和 土壤 
类 型 的 典型 值 进行 土壤 属性 定量 模拟 ,所 得 土壤 
类 型 图 精度 均 优 于 传统 方法 “” ” 。 还 有 学 者 将 
土壤 一 环境 推理 与 目的 性 采样 相 结 合 , 提 高 了 精 
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细 土 壤 普 查 的 效率 和 数字 土壤 制图 的 精 
jag 19.20,22.25-26] 但 由 于 土壤 自身 的 空间 变异 性 及 
其 与 环境 条 件 间 的 高 度 非 线 性 , 现 有 推理 模型 的 
稳定 性 较 低 ,难以 在 不 同 区 域 进行 直接 推广 。 

新 疆 巴 音 布 鲁 克 草原 位 于 天 山南 坡 中 部 ,是 我 
国 最 大 的 荒漠 区 亚 高 山高 寒 草 原 ,重要 的 畜牧 生产 
ERa] 。 区 域内 典型 的 高 寒 草 原 土 壤 发 育 特征 
及 其 环境 要 素 具 有 鲜明 的 地 域 特点 ,而 基于 土壤 一 
环境 的 协 变 关 系 开展 草原 土壤 属性 空间 制图 的 定量 
研究 还 相对 较 少 。 因 此 ,本 研究 在 巴 音 布鲁克 典型 
高 寒 草 原 内 选择 研究 区 ,运用 模糊 聚 类 .GIS 和 目的 
性 采样 方法 ,构建 土壤 一 环境 推理 模型 ,以 土壤 pH 
值 为 例 模 拟 区 域 土壤 属性 的 空间 分 布 , 并 与 OK 和 
IDW 等 传统 方法 进行 比较 。 研 究 将 有 助 于 深入 理 
解 新 疆 高 寒 草原 土壤 发 生 、 演 变 过 程 与 环境 要 素 的 
关系 ,丰富 草地 土壤 学 科 的 案例 研究 ,同时 可 提高 山 
区 土壤 调查 的 效率 和 精度 。 


1 研究 区 概况 


研究 区 位 于 新 疆 巴 音 郭 杨 绽 古 自治 州 和 静 县 巴 
音 布鲁克 镇 境内 ,面积 约 为 4 km ,地 势 北 高 南 低 ， 


有 具有 相对 平缓 的 山 珍 , 较 陡 的 山坡 和 宽 平 的 沟谷 , 海 
IRH 2 603 ~2 995 m 之 间 ( 图 1)。 该 区 域 年 平均 气 
温 -4.8 %C ,1 月 最 低 气 温 可 达 -48 °C ,7 月 最 高 气 
s ,年 降水 量 平均 


温 30.5 °C ,年 平均 风速 2.7 m ， 
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276.2 mm ,年 莹 发 量 高 达 1 022.9 ~1 247.5 mm ,年 
日 照 2 466 ~2 616 h ,热能 562.8 kJ .cm 一，a ,全 
年 积 雪 日 达 150 ~ 180 d ,无 绝对 无 霜 期 , 属 典型 的 高 
寒气 候 ”" 。 研 究 区 土壤 主要 为 亚 高 山 草原 土 ,地 表 
植被 覆盖 度 较 高 ,以 高 寒 草原 类 为 主 ,主要 包括 紫花 
针 茅 (Sipa purpurea) 、 线 叶 苹 草 ( Kobresia capillifo- 
lia) IRF ( Polygonum viviparum ) 等 。 


2 研究 方法 


研究 将 首先 辨识 并 空间 量化 区 域内 土壤 一 环境 
协 变 因子 ,建立 基础 要 素数 据 集 , 采 用 模糊 C 均值 
RX (Fuzzy C-means Clustering , FCM ) 方 法 提取 环境 
组 合 ,在 各 组 合 典型 点 野外 采样 获取 表层 土壤 pH 
值 ,最 终 运用 土壤 一 环境 相似 性 推理 方法 实现 研究 
区 表层 土壤 pH 空间 分 布 制图 ,并 与 传统 空间 插值 
方法 进行 比较 。 
2.1 环境 要 素 集 的 建立 

在 土壤 发 育 过 程 中 ,气候 、 母 质 \ 时 间 、 地 形 等 因 
素 都 可 以 影响 土壤 的 pH 值 。 然 而 在 景观 和 小 流域 
尺度 上 ,地 形 起 到 主要 作用 ” ,其 通过 影响 母 
质 和 水 热 条 件 再 分 配 ” ,进而 导致 pH 值 的 空间 变 
化 。 由 于 本 研究 区 范围 较 小 ,气候 和 母 质 要 素 可 视 
为 均一 ,研究 中 未 予 考虑 。 此 外 ,区 域 地 质 年 龄 的 参 
数 化 难度 较 大 ,而 半 定 量 、 基 于 样 点 的 分 析 又 难以 形 
成 具有 足够 精度 的 栅 格 数据 , 故 研究 中 尚 不 考虑 
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1 研究 区 位 置 及 样 点 分 布 


Fig.1 Location of the study area and field sampling points 
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时 间 因 素 。 因 此 ,本 研究 将 区 域 地 形 特 征 作 为 土壤 
属性 的 主要 协 变 因子 ,参照 土壤 专家 的 意见 ,选取 高 
程 \ 坡 度 \ 坡 向 \ 沿 训 面 曲率 \ 沿 等 高 线 曲率 和 地 形 湿 
度 指 数 6 个 因子 构建 土壤 一 环境 因子 空间 数据 集 
(图 2) 。 基 础 数据 为 分 辨 率 为 30 m 的 DEM, WE 
因子 在 其 基础 上 基于 AreGIS10.0 平台 计算 而 成 。 
2.2 ”基于 模糊 聚 类 的 环境 组 合 获取 

采用 PCM 方法 得 到 的 环境 因子 组 合 代表 了 土 
壤 属 性 在 空间 变化 的 不 同 模式 ,并 以 模糊 隶属 度 的 
形式 表达 ,根据 隶属 度 值 的 高 低 判别 环境 因子 组 合 
中 心 位 置 以 及 组 合 间 的 过 渡 区 ”” 。 该 方法 通过 计 
算 多 属性 空间 中 的 距离 ,使 模糊 聚 类 达到 最 优 , 并 以 
此 确定 各 数据 点 的 隶属 度 “” 。 模 糊 目 标 函 数 如 
F: 


rE > > (ua) (1) 


di, = | yi 2% lla = (ne 0) AC -%) (2) 
式 中 :7 为 模糊 聚 类 的 隶属 度 和 矩阵 ;2 为 聚 类 中 心 
集 ;Y 为 环境 要 素数 据 集 ;n 为 了 中 数据 的 个 数 ;c 
为 聚 类 类 别 的 个 数 ;m 为 加 权 指 数 ;di 为 数据 点 y, 
到 中 心 v; 的 加 权 距 离 ;wi 为 第 个 数据 点 属于 第 i 
类 的 隶属 度 ;4 为 距离 权重 矩阵 ;J 为 模糊 分 类 误 
差 。 


yeY 


FCM 通过 分 割 系数 (Partition coefficient, F ) 和 
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PRELITI ( Normalized entropy , H) 两 个 指标 随 聚 
SCY AE AG TE OUR E RERA BEE 
聚 类 数 的 增加 , 当 分 割 系数 从 (c -1) 到 。 的 递减 量 
比 从 < 到 (c+l) 的 递减 量 小 时 ,可 认为 此 时 的 聚 类 
结果 较 稳定 ,对 应 的 聚 类 数 为 最 佳 '] 。 
2.3 典型 点 土壤 信息 提取 与 土壤 属性 制图 

将 各 环境 因子 组 合 中 隶属 度 最 大 值 所 在 位 置 作 
为 该 组 合 类 的 典型 点 ,构成 土壤 与 环境 的 解释 集 。 
野外 GPS 定位 后 ,在 各 典型 点 采集 土壤 表层 样本 
(0~20 cm) ,经 风干 去 杂 \ 磨 细 \` 过 得 等 处 理 后 , 采 
用 电位 法 测定 土壤 样本 pH 值 ,测定 时 土 水 比 为 
1:2.5, 每 个 土 样 重复 测定 3 次 , 结果 取 平 均 
Ea, 

环境 因子 组 合 类 代表 了 土壤 属性 在 空间 变化 的 
不 同 模式 ,其 模糊 隶属 度 图 可 表达 区 域内 各 像素 对 
其 典型 模式 的 相似 程度 。 基 于 环境 因子 组 合 的 模糊 
隶属 度 和 典型 点 的 土壤 属性 值 ,选择 加 权 平 均 模 型 
推理 研究 区 表层 土壤 pH 值 。 计 算 公式 如 下 . 


V; = > 


式 中 :WV 为 区 域内 未 知 像素 点 (i, J) 的 土壤 属性 值 ; 
VV 为 典型 点 对 应 模式 的 典型 土壤 属性 值 ;5; 为 像 
RU, 门 对 应 环境 组 合 类 的 模糊 隶属 度 值 ;n 为 环境 
组 合 类 总 数 。 


(3) 
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图 2 研究 区 环境 因子 分 布 


Fig.2 Distribution of derived environment parameters of study area 
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2.4 精度 检验 

为 了 检验 上 述 方法 所 得 土壤 pH 空间 分 布 图 的 
精度 ,研究 共 采 集 35 个 独立 样 点 建立 了 验证 集 ( 图 
1) 。 其 中 ,均匀 采样 点 16 个 ,按照 500 m x 500 m 等 
间隔 布置 样 点 ,用 来 检验 区 域 土壤 属性 的 整体 模拟 
精度 ; 横 截 面 采 样 点 10 个 ,设计 为 横 穿 山坡 ,沟谷 的 
线路 , 即 通过 采样 点 在 较 短 的 距离 内 穿越 研究 区 主 
要 的 地 貌 类 型 变化 ,检验 预测 土壤 pH 图 是 否 可 以 
较 好 地 反映 土壤 信息 的 微 地 形 渐变 特征 ;垂直 带 采 
样 9 个 ,以 检验 土壤 pH 值 沿海 拔 梯度 变化 的 模拟 
效果 。 验 证 样本 的 野外 采集 方法 和 实验 室 处 理 分 析 
流程 与 解释 集 相 同 。 选 取 平均 绝对 误差 (Mean Ab- 
solute Error, MAE) 、 均 方 根 误差 (Root Mean Squared 
Error , RMSE ) 和 一 致 性 系数 (Agreement Coefficient, 
AC) 对 推理 结果 进行 定量 评价 ,并 与 OK 和 IDW 方 
法 进行 比较 。 


3 结 R 


3.1 训练 样本 集 构建 

研究 采用 5 种 不 同 的 模糊 权重 对 数据 层 进行 迭 
代 计 算 (m=1.5,1.75,2,2. 25,2.5)'") ,各 权重 下 
均 产 生 聚 类 数 从 2 到 20 的 不 同 聚 类 结果 ,基于 公式 
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(1) A 2) 7 aE RARE n 增 大 产生 的 差 值 。 
结果 显示 , 当 m=2.c=9 Ih} AMIR H. F ke Be 
前 后 两 类 变化 小 且 趋 于 稳定 , 故 本 研究 中 选择 最 佳 
类 别 数 为 9 ,加权 指数 为 2 ,进而 得 到 各 环境 因子 聚 
类 组 合 模糊 隶属 度 分 布 。 图 3 列举 了 9 种 组 合 中 的 
3 组 ,classl 代表 坡 脚 和 沟谷 区 域 , class2 代表 海拔 
相对 较 低 .坡度 较 缓 的 平原 部 位 ,class3 代表 海拔 相 
对 较 高 的 坡 肩 部 位 。 随 后 将 各 环境 因子 组 合 中 隶属 
度 值 最 大 处 的 样 点 作为 土壤 属性 推理 的 解释 样 点 ， 
构成 解释 集 (图 1)。 
3.2 ”表层 土壤 pH 值 推理 制图 

对 典型 样 点 集 进行 一 般 性 统计 ,pH 值 范围 为 
7.287 ~ 8. 263, 平 均值 为 7. 595, 变异 系数 为 
3.32% ,表明 研究 区 土壤 为 典型 的 弱 碱 性 土壤 。 根 
据 土壤 一 环境 的 空间 隶属 度 分 布 ,结合 环境 要 素数 
据 库 和 实验 室 数据 ,由 公式 (3 ) 得 到 研究 区 0 ~ 20 
om 表层 土壤 pH 值 空间 分 布 (图 4a) 。 结 果 显 示 研 
究 区 表层 土壤 整体 呈现 弱 碱 性 ,pH 值 随 地 形 连 续 变 
化 ,变化 范围 为 7.170 ~ 8.186。 与 采用 OK 和 IDW 
的 制图 结果 (图 4b、 图 4c) 相 比 ,FCM 方法 所 预测 的 
土壤 pH 分 布 图 较 好 地 反映 出 表层 土壤 pH 对 地 形 
变化 的 响应 ,显示 出 研究 区 表层 土壤 pH 空间 分 布 
的 更 多 细节 ,对 于 土壤 模拟 结果 更 加 接近 实际 。 
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到 3 环境 组 合 的 模糊 隶属 度 分 布 (以 classl \class2 和 class3 为 例 ) 


Fig.3 Fuzzy membership maps of combinations of environment factors based on FCM( classl , class2 and class3 ) 
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图 4 表层 土壤 pH 值 空间 分 布 模拟 结果 (a-FCM,b - OK ,c - IDW) 
Fig.4 Simulation results of spatial distribution of pH in surface soil(a-FCM,b-OK,c-IDVW) 
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朱 WE ETR R S SEL RE SER pH 值 空间 制图 


3.3 精度 评价 

结果 显示 , 除 个 别 样 点 偏差 相对 较 大 外 ,推理 结 
果 与 实测 结果 总 体 基本 吻合 ,距离 聚 类 中 心 较 近 且 
地 形变 化 较 显 著 的 区 域 ,模拟 精度 相对 较 高 。 全 部 
验证 点 的 整体 平均 绝对 误差 和 均 方 根 误差 分 别 为 


对 较 高 。 这 也 反映 模型 应 用 在 环境 梯度 变化 大 的 地 
方 具有 更 好 的 有 效 性 。 

对 比 土壤 一 坏 境 推理 方法 和 其 它 两 种 方法 的 模 
拟 精度 ,可 以 得 出 :OK 法 和 IDW 法 的 MAE 和 RMSE 
较为 接近 , 均 高 于 FCM 法 ,说 明 本 研究 所 采用 的 土 


0.14 和 0. 19 ,不 同 采样 方式 表现 出 的 模拟 精度 特征 
有 所 差异 。 横 穿山 坡 沟谷 的 线路 的 横 截 面 采 样 点 
的 模拟 值 与 实测 值 最 为 接近 ;垂直 带 采 样 方式 次 之 ， 
均匀 采样 点 的 模拟 精度 相对 最 低 。 横 截面 采样 点 在 
较 短 的 距离 内 穿越 了 研究 区 内 地 貌 类 型 变化 相对 较 
大 的 山坡 沟谷 ,因此 所 建立 的 推测 模型 能 更 好 的 捕 
换 地 形变 化 较 大 处 的 pH 值 变化 ,所 以 模拟 精度 相 


培 一 环境 推理 方法 所 得 预测 结果 精度 高 于 其 它 两 种 
方法 (图 5 和 表 1)。 此 外 ,将 3 种 方法 在 验证 点 处 
的 预测 值 与 实测 值 进行 比较 ,结果 显示 传统 插值 方 
法 预测 值 的 回归 决定 系数 分 别 为 0.22 和 0. 29 , 低 于 
FCM 法 的 0.47, 可 见 FCM 法 的 预测 精度 高 于 OK 插 
值 和 IDW 插值 方法 。 同 时 ,一 致 性 系数 显示 土壤 
环境 推理 方法 的 稳定 性 较 高 。 
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图 5 不 同 模拟 结果 比较 (a-FCM,bp-OK,ec-IDVW) 


Fig.5 Comparison of different simulation results(a-FCM,b-OK,c-IDW) 


表 1 不 同方 法 的 精度 评价 结果 对 比 


Tab.1 Accuracy evaluation of different methods 


且 受 限于 数据 的 空间 量化 可 能 性 与 可 获取 性 ,本 研 
究 中 对 于 少量 野外 样本 的 确定 主要 基于 数字 高 程 模 
型 及 其 衍生 出 的 地 形 因 子 , 未 将 其 它 土 壤 环境 因素 


考虑 在 内 。 如 研究 区 扩大 至 流域 和 区 域 广度 时 , 气 
修 因子 、 母 质 因子 和 人 类 活动 因子 则 是 模拟 时 必须 
要 考虑 的 因素 ,将 有 助 于 样 点 典型 性 模型 模拟 精度 


OK 法 和 IDW 法 是 最 为 常用 的 两 种 空间 插值 方 
法 ,但 它们 只 注重 样 点 的 空间 位 置 和 数据 本 号 ,未 考 
虑 环境 因素 的 作用 。0OLS 法 虽 注 重 环境 因素 的 作 
用 ,忽视 了 空间 位 置 的 影响 ”“: 。 近 年 来 被 广泛 使 用 
的 地 理 权 重 回 归 法 (Geographically Weighted Regres- 
sion, GWR ) ,对 OLS 方法 进行 了 改进 ,同时 考虑 了 


模拟 方法 ”验证 点 类 型 MAE RMSE AC 
FCM 总 体 0.14 0.19 0.81 
均匀 采样 0.33 0.26 
横 截面 采样 0.08 0.09 
垂直 带 采样 0.10 0.13 
OK 总 体 0.23 0.29 0. 68 的 提高 。 
均匀 采样 0.51 0.36 
横 截 面 采样 0. 15 0.16 
垂直 带 采样 0.22 0.25 
IDW 总 体 0.22 0.28 0.72 
均匀 采样 0.51 0.35 
横 截 面 采 样 0.14 0.15 
垂直 带 采 样 0.18 0.23 
4 讨 论 


目的 性 采样 过 程 中 , 样 点 的 选 定 是 土壤 属性 模 
拟 精度 的 重要 前 提 。 景 观 尺度 下 ,决定 土壤 变化 的 
主要 因子 是 地 形 和 水 文 状况 '" ,尤其 在 较 小 空间 范 
围 里 的 精细 土壤 调查 中 ,地 形 要 素 可 代表 土壤 形成 
与 发 展 的 主要 影响 因素 ”。 由 于 研究 区 面积 较 小 ， 


插值 点 的 空间 位 置 和 环境 要 素 的 作用 , 常 被 用 于 对 
土壤 属性 模拟 和 空间 非 稳定 关系 的 探讨 “-” 。 但 
其 本 质 上 是 局 部 地 理 回归 的 结合 , 随 着 被 估计 的 局 
部 模型 越 来 越 多 ,显著 性 水 平 测试 有 意义 的 概率 会 
越 来 越 大 。 此 外 ,不 恰当 的 带宽 参数 会 引起 回归 系 
数 出 现 强 烈 的 空间 变异 或 平滑 ,被 估计 的 回归 系数 
的 空间 模式 被 人 为 地 引入 ,影响 实际 回归 效 
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PS 。 上 述 方法 中 样本 的 数量 对 模拟 精度 影响 
较 大 , 当 样本 较 少 时 ,往往 难以 建立 可 靠 的 模型 。 基 
于 土壤 一 环境 协 变 关系 和 模糊 聚 类 的 建 模 方法 针对 
典型 环境 类 型 进行 目的 性 采样 ,能 够 在 有 限 的 样本 
数量 下 保证 一 定 的 模拟 精度 ,因此 在 兼顾 效率 与 精 
度 时 ,具有 相对 的 优越 性 。 

基于 DEM 的 数字 地 形 分 析 , 作 为 定量 提取 各 种 
地 形 信 息 的 有 效 手段 ,已 被 广泛 应 用 于 水 文学 .地貌 
学 .土壤 学 .生态 学 等 诸多 研究 领域 “1 。DEM 分 辨 
率 的 变化 会 直接 影响 地 形 分 析 的 结果 “1 ,在 实际 
应 用 中 可 通过 对 比 选择 适宜 的 DEM AR, 
研究 所 用 DEM 数据 的 空间 分 辩 率 为 30 m ,提取 的 
各 地 形 因 子 均 在 此 基础 上 完成 。 但 如 采用 不 同 分 辨 
率 的 地 形 数据 ,在 运用 模糊 聚 类 提取 环境 组 合 时 会 
存在 怎样 的 差异 目前 尚 缺 乏 研 究 , 今 后 可 通过 模拟 
精度 比较 分 析 得 出 土壤 属性 推理 的 适宜 分 辨 率 和 最 
佳 推理 空间 。 此 外 ,本 研究 属于 小 尺度 的 土壤 属性 
模拟 , 现 有 研究 中 已 有 不 少 学 者 在 不 同 空间 尺度 开 
展 了 相关 研究 ”3 ,但 土壤 一 环境 推理 的 过 程 和 结 
果 是 否 具 有 尺度 效应 ,如 何 进行 多 尺度 推演 等 ,有 待 
深入 研究 。 


5 结 it 


本 研究 以 新 疆 巴 音 布 鲁 克 山 区 典型 高 寒 草 原 为 
对 象 ,选取 高 程 .坡度 、 坡 向 、 沿 剖面 曲率 、 沿 等 高 线 
曲率 、 地 形 湿 度 指数 等 地 形 因 子 进行 模糊 聚 类 ,获得 
研究 区 影响 土壤 属性 的 典型 环境 因子 组 合 , 结 合 
的 性 采样 ,运用 土壤 一 环境 关系 模拟 人 研究 区 表层 土 
SE pH 值 空间 分 布 ,并 与 传统 空间 插值 方法 进行 比 
较 。 结 论 如 下 : 

(1) 基于 构建 的 基础 数据 集 ,采用 模糊 c 均值 
聚 类 方法 得 到 了 9 类 环境 因子 组 合 , 并 在 各 组 合 隶 
属 度 值 最 高 的 中 心 设 置 了 18 个 典型 点 ,运用 土壤 一 
环境 关系 模型 模拟 得 到 研究 区 表层 土壤 pH 值 空 间 
分 布 图 ,区 域 变 化 范围 7. 170 ~ 8.186 之 间 , 整 体 呈 
弱 碱 性 。 

(2) 模糊 聚 类 方法 通过 较 少 的 典型 样 点 获取 了 
区 域内 土壤 属性 的 连续 性 空间 分 布 , 提 高 了 采样 的 
有 效 性 ,模拟 结果 较 好 地 反映 了 表层 土壤 pH 对 地 
形变 化 的 响应 ,显示 出 研究 区 表层 土壤 pH 空间 分 
布 的 更 多 细节 。 

(3) 经 35 个 独立 样本 检验 ,人 研究 区 土壤 pH {EL 
模拟 结果 与 实测 值 基本 吻合 , 且 模 拟 精度 高 于 普通 


克 里 格 法 (OK) 和 反 距 离 权 重 法 (IDW) ,是 一 种 相 
对 低 成 本 、 高 效率 的 数字 土壤 制图 方法 ,在 草地 土壤 
调查 中 具有 较 好 的 应 用 潜力 。 
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朱 _ 短 等 :基于 模糊 聚 类 的 新 疆 典 型 高 寒 草原 土壤 pH 值 空间 制图 


Predictive soil pH mapping based on fuzzy clustering in 


typical alpine grassland of Xinjiang 


ZHU Lei ， SHENG Jian-dong'”, JIA Hong-tao 
(1 College of Grassland and Environment Sciences , Xinjiang Agricultural University , Urumqi 830052 , Xinjiang , China; 
2 Xinjiang Key Laboratory of Soil and Plant Ecological Processes , Urumqi 830052 , Xinjiang , China) 


Abstract; Obtaining spatial distribution of grassland soil properties accurately and efficiently can provide basic 
information and reference for grassland resource management. Compared with the traditional soil survey method, the 
soil environment inference model based on fuzzy logic can improve the efficiency of field sampling and accuracy of 
predictive mapping, which has been widely used in digital soil mapping. However ,due to the soil’s spatial variation 
and its non-linearity with the environmental conditions, the stability of the existing models are relatively low. The 
models have not been applied in alpine meadow area. In this study ,fuzzy C-means clustering (FCM) was used to 
predict soil pH in the surface layer of grassland soil within a 4 km’ area in Bayanbulak District, Xinjiang Uyghur 
Autonomous Region, China. Six terrain factors, including elevation, slope , aspect, planform curvature , profile curva- 
ture and topographic wetness index , were clustered. Fuzzy membership of 9 groups of environmental factors were de- 
rived to position 18 soil samples in the area with membership larger than 0.9. Then pH distribution was predicted 
with fuzzy membership model. The pH value of study ranged from 7. 170 to 8. 186 and was consistent with the meas- 
ured values. The mapping results reflected continuous changing of soil properties with terrain changing. There were 
35 individual soil samples (16 equal-interval sampling points ,9 cross-sectional sampling points and 10 sampling 
points according to altitude) collected as validation data set. The agreement coefficients between observed values 
and predicted values were high, and the accuracy of FCM model is higher than that of Ordinary Kriging method and 
Inverse Distance Weighted method. FCM and purposive sampling for digital soil mapping is also suitable for small- 
scale region. This approach is an efficient digital soil mapping method with satisfactory prediction precision using 
less samples. It could be possibly applied to the areas with the similar landscape conditions. 


Key words: fuzzy clustering; purposive sampling; pH; Bayanbulak 


